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We consider the neutral, one-dimensional Falicov-Kimball model at zero temperature in the limit 
of a large electron-ion attractive potential, U . By calculating the general n-ion interaction terms 
to leading order in 1/U we argue that the ground-state of the model exhibits the behavior of a 
complete devil's staircase. 

. . PACS numbers: 05.30.-d, 74.25.Dw, 71.30.+h, 

0^ ' In this letter we study the ground-state phase diagram of the one-dimensional Falicov-Kimball model. This model 
was proposed to describe metal-insulator transitions [ and has since been investigated in connection with a variety 
of problems such as binary alloys [ ^ , ordering in mixed- valence systems [ || , and the formation of ionic crystals [ Q . 
O It is the latter language we shall use here, considering a system of static positive ions and mobile spinless electrons. 
2^ , The model comprises no electron-electron or ion-ion interactions but an on-site electron- ion attraction, —U . 
We write the Falicov-Kimball model in the form 

2l : n = t Y,{a]a,+^ + a]+ia,) - U Y,{s,a]a, - 1/2) 
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' where al {ai) denotes the fermionic creation (destruction) operator for a spinless electron, Si is equal to 1 (0) if site i 
C<\ • is (un)occupied by an ion, t is the hopping integral for electrons, /x^ and /Xe are the chemical potentials for ions and 
^-H ' electrons respectively and C/ is a positive constant corresponding to the ion-electron attractive energy. The choice of 
a positive U is not restrictive since the transformation {U —U; fj,i — > —fJ-f, Si ^ 1 — Si} maps the Hamiltonian (|l|) 
onto the same system with U negative. 

The ground state of the system is chosen by minimizing the energy per site over all possible ionic arrangements. 
The structure of the ground states differs significantly depending on whether U is large or small compared to t. In the 
^ ' first case the electrons are essentially localized near the ions and the latter tend to be as far apart as possible while, 
for large t/U, the delocalization of electrons favors the formation of clusters of ions [ ||. In this letter we consider the 
(~| ' case where U is very large compared to other parameters in (|^), and treat t/U as a perturbative parameter. 
O ■ Despite the simplicity of the Falicov-Kimball model the determination of the ground state is far from trivial. 
^ Numerical results [ ^ have suggested that in the neutral system, where the number of electrons and ions are equal, 
^ I a large number of modulated phases appear as ground states. In 1989 Barma and Subrahmanyam studied the phase 
diagram of the model by mapping it onto an Ising system [ 0| . They showed that the phases appearing at the first 
few stages of a perturbative analysis could be described in terms of a simple branching rule, hence suggesting that 
' the complete phase diagram might display a devil's staircase. A different approach to the large-?7 limit was later 
" " ' introduced by Gruber et al. [ ^ ] who considered the model as a set of ions with interactions mediated by the electrons. 
They calculated the two-ion interaction to leading order in t/U on the basis of which they argued that the ion spacing 
is constant in the ground state. 

Here we show that a full determination of the ground state requires a calculation of the general m-ion interactions. 
These are obtained to leading order in t/U using Green's function techniques. Then, using arguments first introduced 
by Fisher and Szpilka [ ||, we deduce the existence of a devil's staircase in the neutral Falicov-Kimball model. 

The phase diagram for t = is shown in Fig. |l|. All the phase boundaries in the figure are multidegenerate in 
that any phase obtained by mixing the two neighboring phases is degenerate on the boundary. Our aim is to study 
systematically how this multidegeneracy is lifted as t/U increases from zero. 
It is convenient to introduce the variables 

h=(p, + ^l,)/2 , (2) 

A = (m. - Me)/2 . (3) 

U is assumed to be much larger than any physical parameter in (|l|) and therefore A/U ^ 1. This restriction on A has 
the important consequence of fixing the total number of electrons equal to the total number of ions, and throughout 
the rest of the paper, we will implicitly consider a neutral system, J^i — J2i where Ui — a|ai. 
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When moving along the line /Xg = fii in Fig. yjone notices that, for negative values of /i, the ground state corresponds 
to an empty lattice (n^ = Si = 0). On the other hand, for h positive Ui = Si = 1. The point h = lies on the multi- 
degenerate phase boundary where all phases associated with an arbitrary spacing of the ions are degenerate. To 
distinguish between the different degenerate states it is convenient to introduce the labelling (ni, n2, ..rim) to denote a 
phase consisting of ions whose separations (measured in lattice spacings) repeat periodically the sequence ni, n2, .., rim- 
(Hence the phases rii = Si ^ 1 and rt,; = = can be described as (1) and (oo) respectively.) 

The multidegeneracy encountered on the phase boundaries of Fig. |l| is due to the absence of interaction between 
the confined electrons. It is natural to expect that, for t/U ^ 0, the hopping of electrons will introduce an effective 
coupling between the ions, thus providing a mechanism for the removal of the degeneracy. This intuitive picture can 
be formalized using the defect-defect interactions introduced by Fisher and Szpilka [ |^. In the present context, a 
defect corresponds to an ion. Following [ ^ the energy per lattice site of phase (ni, ^2, n^) can be written as 
£^{ni,...,n„) = ^totlY^=\ "■i; where 

rn m 

Etot = mo ^^V2(ni) ^^Vi,{ni,ni^x) ^ - (4) 

where a is the creation energy of an isolated ion, V2(a;) denotes the effective interaction between two ions at a distance 
a;, V'iix^y) the interaction of three ions with spacings a;, y, and so on. Although, for simplicity, we refer to the ion 
creation energy and ion-ion interactions, it must be borne in mind that each ion is associated with an electron. 

When t = the electrons are confined to the ions. In this case a is readily shown to be equal to —2h. For small tjU 
we expect each electron to be localized in a region around the associated ion. Using standard perturbation theory 
one can obtain a to leading order va^tjU 

a^-2h-2t^/U + 0{t'^/U^) . (5) 

The leading order corrections to a are associated with virtual a process in which the electron hops to the site to the 
immediate right (or left) of the ion and back again. 

The general ion-ion interaction term, T^(ni,n2, ...,nm_i), can be obtained, at least in principle, through a recon- 
nection formula [ ^. In terms of the four different ionic configurations shown in Fig. ^ this formula is 

Vm{ni,n2, ...,n„_i) ^ Ea - Eb - Ec + Ed ■ (6) 

In the absence of electron hopping, Eq. (^ gives Vm = for all values of m. Our aim is to calculate Vm to leading 
order in t/U . A simple way to obtain the leading-order contribution to the interaction between two ions occupying 
sites and n relies on perturbation diagrams in which the matrix element of talaj is represented by an arrow from j 
to i and the ordering of matrix elements follows the height on the page. First note that diagrams involving disjoint 
sets of sites do not contribute to the ground state energy. In addition, diagrams involving fewer than 2n electron hops 
will not contribute to V2{n) because, as illustrated in Fig. the contribution of every such diagram in configuration 
A will be cancelled by a counter-diagram in configuration B (or C). 

It is then apparent that the leading-order contribution to V2(n) is due to diagrams where the sites and n are 
just connected by 2n hoppings, as in Fig. ||b, and is proportional to The proportionality factor can be 

calculated by summing the contributions from all relevant diagrams 

V2in) = 2nt2"/;72"-i + C'(<2"+2/[/2n+i) . (7) 

Fisher and Szpilka [||] showed that, for systems where the ion-ion interactions, Vm, decay sufficiently rapidly with 
the defect spacings, a knowledge of the sign and convexity of V2(n) can provide a considerable amount of qualitative 
information about the phase diagram of the system. Their analysis can be applied in this context since the Vm decay 
exponentially with the spacings of the two outermost ions (because the ion-ion interaction is mediated by the nearest- 
neighbor hoppings of electrons). Therefore, as a first approximation, we shall analyse the phase diagram neglecting 
interactions that involve more than two ions. Higher-order interactions will then be included successively to resolve 
the finer details of the phase structure. 

In the two-ion interaction approximation the ground-state configurations correspond to equispaced electron-ion 
pairs. Since V2(n) is always positive and convex, as h is varied from positive to negative, n increases monotonically 
in steps of one lattice spacing [ |^ , giving rise to the infinite sequence of phases 

(1) ^ (2) ^ ... ^ (c^) . (8) 

The phase (n) is stable over a region of width 
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The original multidegeneracy is not completely lifted by V2{n) because, on the boundary between two phases, (n) 
and (n + 1), all mixed phases where the ions can be separated by distances n or n + 1 are still degenerate. To determine 
the finer structure of the phase diagram it is necessary to consider the effect of higher-order ion interactions. These are 
not easily obtained using the simple method outlined above since, when there are more than two ions, it is extremely 
difhcult to keep track of the energy denominators associated with the different orderings of the electron hoppings. 
However, this problem can be circumvented by using Green's function techniques. We first illustrate how this method 
can be used to reproduce the result for V2{n). 

To calculate in Eq. (^) consider a system of n + 1 sites with ions at sites and n. The single-particle energies 
are determined by the eigenvalues of the x(n -I- l)-dimensional matrix, M, where 



(10) 



and the other matrix elements are zero. Two of these energies occur near — [/, and these are the ones we want to sum 
over. So we write 
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where the contour F encloses the region near z = —U and T is the identity matrix. To evaluate the trace we expand 



the matrix inverse in Eq. ( |ll| ) in powers of the fs. Define a perturbation Vij = t{6j^i^i + 6j^i-i). Then 
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where Ga — [z — Ma]^^. Terms which are odd order in t cannot contribute to the trace. 

In this expansion one sees that, if i is not an end site, in order to involve all the fs the matrix elements must start 
at i, say, then increase to the highest number site (n), then decrease to the lowest number site (0) and finally increase 
back to the original value i. Alternatively, the matrix elements could initially decrease. If i = or n however, note 
that the index can only initially increase or decrease respectively. So to leading order 



zJ —M.\ ~ GiGooGiiGn,n Y\_ ^'ii Y\. ' 



(13) 



where = 1 if i = or i = ri and Ci — 2 otherwise. The product over G"s does not include the end sites, because 
these, in general, only appear once. The starting site appears an extra time and gives rise to the prefactor Ga. The 
term of order t^" in Eq. ( 1^) is 



Tr(0X-M)~^ « t 



2n 



2n-2 



{z + Ufz^'^-^ (z + U)^z^^-^ 
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Here the first term includes Ci and C„+i, both of which are unity. The factor 2n — 2 comes from Y^^=2 ^i- Substituting 
( p^ ) in (O) and calculating the integral using residues gives (here and below we give the expressions only to leading 
order in t/U) 



2m ./r 



2n 



-I- 



2n- 2 
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(2n-2)i2"/[/2"-i . 
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Next, to use the reconnection formula (|^), we need to repeat the same calculation when one of the end ions is removed 
(corresponding to configurations B and C in Fig. In this case 
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Thus the perturbative contributions to V2, denoted Eb and Ec, are Eb = Ec = —t^^/U'^"'~^. Note than when both 
ions are removed there are no longer any levels near —U. Hence Eu = and use of the reconnection formula (|^) gives 
V2(n) = 2nt^"/[/^"~^, in agreement with the expression (^. 

The method outlined above can be extended to calculate the m-ion interaction Vm for m > 2. As we shall show 
below, Vm{ni,n2, ...nm-i) depends, to leading order, only on the separation of the two outermost ions in configuration 
A, n — The result is 

V (n) ^ (16) 

" (2m - 3)!(2n - 2m + 3)! C/2"-i " ^ 

To prove this consider Eq. (|l3|). Note that m of the diagonal elements of {zl — M)^^ are (z + U)^^; the rest are 
z^^ . If the initial i corresponds to an ion, then a factor (2 + [/)2ni-i appears in the trace; otherwise the factor is 
[z + C/)2™~2_ lY^Q gj.g^ Qase there are m choices for i; two at the end with Ci = 1 and m — 2 in the interior with 
a ^2. Thus 



Tr 



(zJ-M)-i 



^^2n] (2m -2) 



I {z + [/)2m-1^2n-2m+2 

2n - 2m + 2 1 

(z + [/)2m-2^2n-2m+3 | ' 



(17) 



Again we stress that the dependence of ( p^ ) on the position of the m ions in the chain is only through n, the distance 
between the two end defects. Substituting in (ni|) gives 



(2»"2)! 

^ (2m-3)!(2n-2m+l)!C/2»-i ■ ^ ' 

Similarly 

= i^c = ^A(2m - 3)/(2n- 2m + 2) (19) 
Ed = EB{2m~4:)/{2n-2m + 3) . (20) 

Finally the use of the reconnection formula (||) gives for the m-ion effective interaction, Vm , the result (|l^) . It should 
be pointed out that, in principle, the leading order expression ( p^ ) could be dominated by neglected terms of higher 
order in t/U if n is sufficiently large (for fixed t/U) However, Gruber et a/. have shown that, for m = 2, 

higher-order corrections to V2 (n) are dominated uniformly in n by the expression (n6|) , provided that t is replaced by 
i = UIVW^ + it"^ — U]/2t. It seems plausible to expect that, upon renormalizing t in ([l^), their conclusion can also 
be extended to m > 2. 

We now consider how higher-order ion interactions modify the phase diagram obtained in the two-ion interaction 
approximation. Consider first V3. This has the effect of partially removing the multidegeneracy on the {n)\{n + 1) 
boundaries by stabilizing the mixed phases {n,n+ 1). This happens because the energy difference 



(2n + l)£^(„,„+i) - n£'(„) - {n + l)£^(„+i) = 

V3{n,n + 1) + Vsin + l,n)- ¥3(71,71) - l/3(n + l,n+ 1) 



is negative. The mixed phase (n, n -I- 1) has an ion density, 2/ (2n + 1), intermediate between the pure phases (n) and 
(n + 1). 

The stability of the two new boundaries appearing at this stage of approximation, namely (n)|(n,n -I- 1) and 
(n,n -|- l)|(n -|- 1) can be determined similarly by considering four-ion interaction terms. Again they are unstable 
to the appearance of the mixed phases {71,71,71 + 1) and (n, n + l,n + 1) respectively. Indeed, since all interaction 
potentials are positive and decay exponentially with the separation of the outmost ions, we can conclude that, at 
every stage of the construction of the phase diagram, the introduction of neglected higher-order interactions will lead 
to the stabilization of mixed phases of increasingly long period. 

To summarize: we have calculated the general m-ion interaction potentials in the neutral Falicov-Kimball model 
to leading order in t/U at zero temperature. We thereby iteratively construct the ground-state phase diagram and 
conclude that the ion density versus chemical potential, h, has the form of a complete devil's staircase. 

Extending the strategy for the iterative construction of the phase diagram to more than one dimension is not trivial 
and is the focus of an ongoing investigation. 
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FIG. 1. The phase diagram of the Falicov-Kimball model for t = 0. 
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FIG. 2. Ionic configurations needed to calculate the m-ion interaction Vm{ni,n2, ...rim-i). In A there are m ions with 
successive separations ni,n2, ...,nm-i- In B the left-most ion is removed; in C the right-most ion is removed; and in D both 
the left-most and right-most ions are missing. 
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FIG. 3. Examples of diagrams that, when the reconnection formula is used, give an (a) zero (b) leading-order contribution 
to the two- ion interaction V2(4). A full circle represents an ion and an arrow denotes the hopping of an electron. The contribution 
to the energy can depend on the order of the arrows which determines the energy denominators, pU , that arise in perturbation 
theory, where p is the number of electrons away from their ions. For example, in configuration A\, six denominators are U and 
one is 2U whereas, in A2, five denominators are U and two are 2U . 
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